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ABSTRACT 

We study the role of spin-down in driving quark deconfinement in the high density core of isolated 
neutron stars. Assuming spin-down to be solely due to magnetic braking, we obtain typical timescales 
to quark deconfinement for neutron stars that are born with Keplerian frequencies. Employing differ- 
ent equations of state (EOS), we determine the minimum and maximum neutron star masses that will 
allow for deconfinement via spin-down only. We find that the time to reach deconfinement is strongly 
dependent on the magnetic field and that this time is least for EOS that support the largest minimum 
mass at zero spin, unless rotational effects on stellar structure are large. For a fiducial critical density 
of 5po for the transition to the quark phase (po — 2.5 x 10 14 g/cm 3 is the saturation density of nuclear 
matter), we find that neutron stars lighter than 1.5M cannot reach a deconfined phase. Depending 
on the EOS, neutron stars of more than 1.5Mq can enter a quark phase only if they are spinning 
faster than about 3 milliseconds as observed now, whereas larger spin periods imply that they are 
either already quark stars or will never become one. 
Subject headings: neutron stars, spin-down, quark deconfinement 



1. INTRODUCTION 

Neutron stars, formed in the aftermath of core collapse 
supcrnovae, are born rotating with spin periods that can 
be down in the millisecond range. Coupled to this fast ro- 
tation, magnetic flux conservation in the imploding pro- 
genitor provides the nascent neutron star with magnetic 
fields in the range 10 1 2 -10 15 G, as confirmed by a var iety 
of pulsar observations ijLvne fc Graham-Smith 11997(1 . A 
customary first approximation is to model the pulsar as 
a rotating magnetic dipole which loses angular momen- 
tum through electromagnetic radiation. Consequently, 
the star spins down and the decreased centrifugal force 
causes the core density to rise. As this occurs, bary- 
onic matter can undergo transitions to exotic phases of 
strongly interacting matter, finally giving way to a phase 
in which quarks are deconfined. Such a deconfined phase 
might explain the most energetic astrophysical phenom- 
ena such as gam ma-ray bursts and cosmic rays from 
compact objects (lOuyed. Keranen fc Maalampil [2005a; 
lOuved. Raon fc VogtJ l2005bll . A phase transition to 
quark matter can result in the formatio n of a quark s tar, 
if certain stability criteria are fulfilled JWeber 112005^ . 

Our main objective in this letter is to obtain accurate 
timescales for spin-down of a neutron star to the decon- 
finement density. Ultimately, we would like to know the 
likelihood of a neutron star to undergo a phase conver- 
sion to a quark star. To this end, we explore 4 differ- 
ent EOS to determine plausible bounds on masses and 
spin-periods that are consistent with a phase transition 
to quark matter in the star's spin-down lifetime. As- 
pects of spin-down by magnetic braking and subsequent 
deconfinement are outlined in section [21 along with the 
approximations made in this work. In section[21 we elab- 
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orate on our computations and the EOS used (see Fig.^l, 
connecting them to results displayed in Tableland Fig- 
ure [2J Section 0] gathers our principal conclusions based 
upon an interpretation of results in Tableland Figures[21 

2. SPIN-DOWN AND DECONFINEMENT 

It is well-known that the energetics of pulsed emission 
as well as unpulsed nebular radiation from a neutron 
star can be generally un derstood as deriving from it s 
rotational kinetic energy [(Manchester fc Tavlor 111977(1 . 
While the rotational power can be calculated from mea- 
surements of the spin-period (P) and its derivative (P) in 
a model-independent way, theoretical calculations link- 
ing surface magnetic fields with the spin-down rate (loss 
of rotational kinetic energy) require a model for the mag- 
netic field configuration. In this paper, we assume a dipo- 
lar time-independent magnetic field for the sake of sim- 
plicity. Although more complicated configurations with 
toroidal components are likely in reality, these will not 
affect our quantitative results which depend more on the 
absolute strength of the magnetic field rather than its de- 
tailed shape. Further, surface magnetic fields are likely to 
decay over a million years or so, and we do not consider 
this complication here since its qualitative impact can 
be easily inferred from our results. We further assume 
that gravitational radiation is a negligible contributor to 
spin down, i.e. our neutron stars are born axisymmetric 
and remain as such, so that the spin-down rate of a neu- 
tron star is entire ly due to dipolar magnetic field braking, 
as gi ven by (e.g. iDeutschl 119551 : iManchester fc Tavlor I 
[HZ3): 
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where B is the equatorial surface magnetic field, / is the 
moment of inertia, R the equatorial radius of the star, 
£1 = 2ir/P its instantaneous angular velocity and c the 
speed of light. In equation JTJ above, the inclination 
angle of the dipole (x) has been averaged over a sphere 
representing a population of fixed mass and magnetic 
field ((sin 2 x) = 2/3). With all quantities expressed in 
cgs units, equation^ can be integrated to yield the time 
for a neutron star to spin down to a particular 17. The 
additional (implicit) input that is required to solve for 
the spin-down time is the EOS, which determines the R 
and / for a given il. 

Computing the spin-down time to dcconfincmcnt 
requires precise knowledge of the details of a phase 
transition to quark matter. Although yet to be de- 
termined exactly from first principles, it is expected 
to occur on general grounds at high baryonic den- 
sity, whose value depends on the E OS employed 
llHeiselberg fc Pandharinande I 120001 iGlendenning I 
2000). If the interface tension between the quark 
and nuclear phases is too high, a first order phase 
transition to a pure quark matter core at very high 
density (~ 10po) is likely, by which time the star 
may become gravitationally unstable and collapse to 
a black hole. Several other phases, such as meson- 
condensed phases or structured nuclear matter can 
intervene at lower densities, softening the EOS. The 
concomitant rapid contraction of the softer matter 
can lead to non-monotonic behavior of the angular 
velocity £7 and cusps in the brakin g index n = ilti/Cl 2 
ijHeiselberg fc Hiorth-Jensen"lll998|) . It is also possible 
that the transition from nuclear to quark matter at high 
density is second-order if the interface tension between 
these phases is small enough. The transition may then 
occur smoothly over a range of densities from (4 — 8)po 
with an incre asing proportion of quark matter in the 
mixed phase ijGlendenning I Il992|) . and a continuous 
variation of the pressure and angular velocities. Our 
working assumptions will be in line with the latter, 
and we choose a threshold deconfinement density of 
5,oo- Obviously, our quantitative results depend on 
this assumed value, but will remain qualitatively true 
irrespective of the transition density one may choose. 
When the central baryonic density reaches the critical 
density of p cr it = 5po> deconfined quark matter will 
appear in the inner regions of t he neutron star. If the 
strange quark matter hypothesis ijltoh 11197(11) holds true, 
it is possible that the entire star converts to a quark 
star in a violent phenom e non termed the Quark Nova 
llOuyed. Dev fc Devll2002l iKeranen. Jaikumar fc Ouvedl 
2004). Aspects and consequ e nces o f the Quark Nova 
are dis cussed in lOuved et al.1 l|2005aj) and lOuved et al.1 



3. NEUTRON STAR MODELS AND EQUATIONS OF STATE 

In order to explore the increase in central density due 
to spin-down, we create sequences of neutron star mod- 
els using the RNS code l|Stergioulas fc Friedman iTlQQfi^ 
developed specifically to treat rapidly rotating neutron 
stars. Models in a particular sequence have the same 
EOS, constant baryon number, increasing central den- 
sity and decreasing angular velocity. We assume that 
the neutron star is born axisymmetric and initially ro- 
tates at Keplerian frequencies. For each sequence, the 
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Fig. 1.— The EOS used (see §3). 

RNS code outputs (among other physical quantities) /, 
tt, and R, which feed into equation^. We emphasize 
that since the magnetic field does not appear explicitly in 
the calculations performed by the RNS code, it does not 
affect the structural evolution of the star or the central 
density at which the transition to quark matter occurs. 
In principle, extremely high magnetic fields exceeding 
10 15 G can affect neutron star structure ijBocauet et al. I 
1995). In our simplified approach, as per equation^), 
the magnetic field only determines the time to reach the 
deconfinement density. 

We employ 4 different EOS which are, in in- 
creas ing order of stiffness: BB B 2 flBaldo et al. I 
Il997lk C (IBethe fc Johnson I 119m APR 
Akmal^jF^andha^irjande&JR.avcnhal] 119981 ) and OBJ 



Jaikumar fc Ouvedll200fil) . The pressure versus density 
curve for each of the above EOS is plotted in Figure Q 
for densities ranging from po to 5po- The OBJ EOS is a 
lot stiffer than the rest at high densities, due to strong 
re pulsion from vector m esons at short distances (details 
in l Jaikumar fc Ouved ikOOfjj) . 

The minimum baryonic mass that is required for a 
given EOS to support spin-down to deconfinement is 
found by the RNS code as the sequence that spins down 
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Fig. 2. — Frequency vs central density for EOS used. All se- 
quences are with constant baryonic mass and for minimum mass 
configurations (see Table . The minimum mass configurations 
are those that can reach a critical density of 5po due to spin down. 
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to p cr it at zero spin. The gravitational mass of such a 
model at zero spin is what we call minimum mass. Fig- 
ure |21 displays the spin frequency of neutron stars with 
minimum mass, as a function of its central density. The 
curves shown are for a 1.46M Q star with EOS BBB2, 
a 1.53M Q star with EOS C, a 1.78M star with EOS 
APR, and a 2.8M star with EOS OBJ (all masses are 
gravitational masses at zero spin) . 

Equation implies that it takes infinitely long to 
spin-down to zero frequency. In practice, we therefore 
compute the time to reach within 1% of p C Ht, as illus- 
trated in Figure For a minimum mass star, we take 
this time to be the maximum time for a star with a given 
EOS to reach deconfinement density. Neutron stars with 
mass lower than the minimum mass will not reach de- 
confinement density due to spin-down, and stars with 
higher masses will reach the critical density sooner, be- 
fore they have spun down completely. Conversely, a max- 
imum mass can also be found, such that any mass beyond 
this would already have central densities exceeding p cr n 
at the Kepler frequency so that it is already in a decon- 
fined state and spin-down causes no further change. 

The minimum and maximum masses as described 
above are listed in Tabled for each EOS. Note that the 
OBJ EOS, on account of its stiffness, predicts a large 
minimum mass in order to spin-down to deconfinement. 
The minimum density at Keplerian frequency required 
to reach deconfinement density (p m i n k) is lowest for the 
OBJ EOS and highest for the APR EOS (these are the 
two stiffest EOS). This is because the moment of inertia 
is so much larger for the OBJ EOS as compared to the 
others, that it more than compensates for the increase in 
angular momentum with increasing stiffness of the EOS. 
This also explains the lower spin frequency at a given 
central density for the OBJ EOS. 

4. RESULTS AND CONCLUSIONS 

The maximum time to reach deconfinement via spin- 
down can be read off from Figure [3] for each EOS for 
magnetic fields in the range 10 12 — 10 15 G. For any of the 
EOS, it takes up to a few hundred years for the star to 
spin down to 0.99p cl n for a 10 12 G magnetic field. This 
time is drastically reduced to about an hour for a 10 15 G 
magnetic field. It is important to emphasize that the 
curves in Figure [3] correspond to the minimum mass star 
in each case, i.e. it takes an infinite time for the magnetic 
field to spin such a star down to /9 C rit- In Table we 
therefore list the time it takes the star to increase its 
central density to within 1% of the critical density (i.e 
0.99/9 cr ;t) due to spin down. 

Figure 0] shows the period as a function of time for 

TABLE 1 

Columnwise: minimum gravitational mass at zero spin, 
maximum gravitational mass at kepler rotation and 
minimum density at kepler rotation that can reach the 

CRITICAL DENSITY (p cr it = 5po) DUE TO SPINDOWN. 



EOS Af min (M ) M max (M Q ) 

Pmin,K (po) 

BBB2 1.46 1.79 3.92 

APR 1.78 2.16 4.12 

C 1.53 1.83 3.88 

OBJ 2.8 3.36 2.52 



each EOS and for different magnetic field strengths. 
Just as before, this is for minimum mass stars and we 
plot the period until the central density is within 1% of 
p cr it. It is noteworthy that the central density does not 
change much after the rotation period slows to about 
-Pmax = 3ms (about 6 ms for the OBJ EOS). This is inde- 
pendent of the magnetic field since that only determines 
the time to deconfinement, not the evolution of period 
with time. From the above limiting value of the period, 
and the minimum mass results tabulated in Tabic ^ for 
our fiducial critical density of 5po, we conclude that only 
stars with mass greater than about 1.5M0 that have a 
spin period P < P max can reach a deconfined phase. 
Conversely, if the star has P > P max as observed now, it 
cannot change its central density sufficiently to undergo 
deconfinement due to spin down, and must either remain 
a neutron star or have been born as a quark star. 

Although the critical density and consequently P m ax is 
not known precisely, the range of currently observed iso- 
lated neutron star masses in conjunction with our study 
implies that such stars can spin-down to a deconfined 
phase only if the deconfinement density is not too much 
above 5po- If the deconfinement density is much larger 
than t his, neutron stars in the currently observed mass 
range llStairs 112001 cannot increase their central densi- 
ties enough via spin-down to reach the threshold, and 
will remain neutron stars for their lifetimes. This con- 
clusion can only be avoided by assuming that all neutron 
stars are actually quark stars and are born as such in the 
aftermath of supernova explosions. Since most realis- 
tic quark matter EOS canno t support a maxim um mass 
greater than 1.6M Q (see e.g. lDev et al. Ill998j) . it is un- 
likely that this is true. We therefore believe that the 
existence of quark stars is predicated upon having a low 
deconfinement threshold (< 5po). Once this threshold 
is pinned down via theoretical studies addressing Quan- 
tum Chromodynamics at high baryon density, we will be 
able to make more definitive statements about the exis- 
tence of quark stars and related phenomena. Until then, 
the possibility that some isolated neutron stars are really 
quark stars cannot be dismissed. 
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TABLE 2 

Time to spin down from Kepler, frequency to the 

FREQUENCY CORRESPONDING TO 1% OF p cr it FOR A STAR WITH 
MINIMUM MASS FOR DIFFERENT MAGNETIC FIELD STRENGTH. 



Time 

EOS M min (M e ) 10 15 G 10 14 G 10 13 G 10 12 G 

BBB2 1.46 96 min 6.7 days 667 days 183 years 

APR 1.78 73 min 5.1 days 505 days 152 years 

C 1.53 89 min 6.2 days 619 days 170 years 

OBJ 2.8 348 min 24 days 2418 days 663 years 
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Fig. 3. — Central density versus time for each of the EOS used with p cr it = 5po- All of these models has the minimum gravitational 
mass (see Table HI , The magnetic field is 10 X G with X = 12, 13, 14, 15. The curves terminate when p is within 1% of p cr it- 
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Fig. 4. — Period vs time for the EOS studied in this paper with p cr it = 5po, with the minimum gravitational masses listed in TablcHland 
for different magnetic fields. The curves terminate when p is within 1% of p cr it- Note that for a given EOS all curves reach the same final 
period since the magnetic field determines only the time to reach deconfinement, and does not enter the hydrostatic equations (see text). 

no. W-31-109-ENG-38. 
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